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Virtual hadronic contributions to the Bhabha process at the NNLO level are discussed. They are substantial for 
predictions with per mil accuracy. The studies of heavy fermion and hadron corrections complete the calculation 
of Bhabha virtual effects at the NNLO level. 



1. Introduction 

Since Loops and Legs 2006 [l|2j . consider- 
able progress in the determination of the virtual 
NNLO corrections to the massive QED Bhabha 
process has been made |l|. 

The Nf = 2 two-loop corrections with heavy- 
fermion insertions have been computed in the 
limit TOg << m? << s,i: with a direct Feynman 
diagram calculation in [T5], and using a factor- 
ization formula that relates massless and massive 
amplitudes in [16]. 

Hadronic contributions have been recently cal- 
culated in [17j . In addition, heavy-fermion cor- 
rections, beyond the nri << s,t limit, have been 
made available in [17118119] : see also [50]. 

Interestingly, the original expectations on the 
necessity of a complete, direct two-loop massive 
Feynman diagram evaluation have not been ful- 
filled yet. After the analytical evaluation of a 
massive planar and a massive non-planar double- 
box diagram (both with seven propagators) in 
[2T] and in [22], there was hope to evaluate all 
the remaining diagrams soon [23] . So far only 
part of the photonic master integrals, namely the 
planar ones, have been evaluated in [53] in the 
limit of small electron mass. Although this limit 
is by far sufficient for experiments, there is still 
space for theoretical developments. 

Finally, it should be stressed that results for 
three gauge-invariant classes of NNLO Feynman 
diagrams have been determined by at least two 
independent groups, relying on different methods: 



iPor earlier literature, see I3I4I5I6I7I8I9I10I11I12I13IT41 . 



• Photonic corrections: computed in [12113] 
and recalculated in [T6l: 



Electron Nf = 1 corrections: computed 
in [T^ and cross-checked in [IS] (with full 
me dependence) and in [16] (small electron 
mass limit); 



• Heavy-fermion Nf = 2 contributions: de- 
termined with two independent methods in 
the limit nii << s,t in [15116] and for any 
mass ruf in [25117] (dispersive approach) 
and in [18|19] (analytical result). 

As far as the full electroweak model is con- 
cerned, the NLO corrections are well known 
[26127] . At NNLO, the terms enhanced by the 
heavy top quark are taken into account [28129] . 
The full dependence of weak corrections on the 
top quark and Higgs boson masses is known 
[30131] and implemented for fermion pair produc- 
tion in e"'"e~-annihilation [32], but is not yet im- 
plemented for Bhabha scattering. Here, we report 
on another class of NNLO corrections, which have 
recently determined: 

• Virtual hadronic NNLO contributions, in- 
cluding both reducible self-energy insertions 
and irreducible vertex and box corrections 

m- 

For a review on the status of Monte Carlo stud- 
ies for Bhabha scattering we refer to [33] ■ 



2. Hadronic Virtual Corrections 

Higher order hadronic corrections to the 
Bhabha scattering cross section can be obtained 
inserting the renormalized irreducible photon vac- 
uum polarization function, 11, in the appropriate 
virtual-photon propagator (see Fig.[T]), 
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where q is the momentum carried by the vir- 
tual photon and S —^ 0+. The vacuum polariza- 
tion function 11 can be represented by the once- 
subtracted dispersion integral [51] . 
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where the production threshold for the interme- 
diate state in 11 is located at q^ = AM"^. We leave 
as understood the subtraction at q^ ^ for the 
renormalized photon self-energy. 

Light-quark contributions get modified by low- 
energy strong- interaction effects, which cannot 
be computed using perturbative QCD. However, 
these contributions can be evaluated using the op- 
tical theorem [35], and relating ImHhad to the 
hadronic cross-section ratio i?had IS] , 

ImHhad (^:) = --i?had(^), (3) 
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Ini Hhad can be computed employing the ex- 
perimental data for i?had iu the low-energy re- 




Figure 1. Two-loop irreducible vertex and box 
topologies for Bhabha scattering with vacuum po- 
larization insertions. 



gion and around hadronic resonances, and the 
perturbative-QCD prediction in the remaining re- 
gions. 

Here we give numerical results through tables 
which have been included also in the first version 
of [360. 

A shape of the Rhad parameterization we use 
is given in Fig. [2l In a forthcoming publication 
[35], containing a detailed description of [T7], we 
will employ a more updated parameterization of 
Rhad [39 40 41j . We just mention that the final 
numbers get modified only slightly and do not 
change qualitatively the situation. 
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Figure 2. Rhadis), based on the parameterization 
given in 



The lower integration boundary in Eq. ([2]) is 
given by M = m^^, where m,^ is the pion mass. 
For self-energy corrections to Bhabha scattering 
at one-loop order this method was first employed 
in [131. 

Finally, we note that contributions to H arising 
from leptons and the top quark can be computed 
directly in perturbation theory, setting M = ruf 
in Eq. ([2]) , where ruf is the mass of the fermion 



^In order to make |17) less technical, we have replaced 
in the final version tables by figures. This means that the 
numbers shown here correspond to the numbers which can 
be read out from I17h see also | 37| . 



appearing in the loop, and inserting the imagi- 
nary part of the analytic result for 11 . 

In the following, we will not consider hadronic 
effects to the running of the coupling constant; 
details can be found in [38] , 

3. Vertex Contributions 

Hadronic irreducible vertex corrections are ob- 
tained through the interference of the vertex dia- 
grams of Figure [T] with the tree-level amplitude. 
Their contribution to the 0{a^) differential cross 
section is given by: 
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where we define vi{x,y) = x"^ + 2y'^ + 2xy and 
V2(x,y) = {x + yY). Here V2 summarizes all two- 
loop fermionic corrections to the QED Dirac form 
factor, whose computation can be traced back to 
the seminal work of [44I45J . The full result can 
be organized as 
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where V2e denotes the electron- loop component, 
see [46] , 

Heavy-fermion and hadronic contributions, in- 
stead, can be evaluated as in Ref. ^47j through 
the dispersion integral 
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where R is given by 
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and where Cf and Qf denote the color factor and 
the electric charge. The two-loop irreducible ver- 
tex kernel function Ky, in the limit of a vanishing 
electron mass, reads as: 
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Here Li2 is the usual dilogarithm and ^2 — 
Li2(l)=^V6. 

In Tables [1] and [2] we show numerical values 
for the various components of V2 (x) of Eq. ^ for 
space-like and time-like values of x {t and s chan- 
nel). For each fermion flavour, we show the result 
obtained through the dispersion-based approach 
(first line) and the one coming from the analytical 
expansion (second line). 

We can see that the latter numbers approach 
the former ones in regions where the analytical 
expansions are expected to become good approx- 
imations. When rai > s, the entry is suppressed. 



4. Box Contributions 

Notice that, unlike the vertex kernel, the irre- 
ducible box kernels are infrared divergent, but, 
analogously to the one- loop box diagrams, they 
have no singularity in the electron mass q In 
order to construct an infrared-finite quantity, we 
combine: (i) Born diagrams interfering with two- 
loop box diagrams and reducible vertices (first 
row in Fig. [3]); (ii) diagrams with a one-loop vac- 
uum polarization insertion interfering with one- 
loop boxes and vertices (second row in Fig. [3]); 
(iii) real single-photon emission diagrams with a 
one-loop vacuum polarization insertion (third row 
in Fig. [3]). The infrared-safe Nf = 2 irreducible 
vertices (see Fig.[T]) and pure self-energy diagrams 
are not included here. 



•^In |38j appropriate simple arguments based on counting 
of logs will be given. 
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Table 1 

Contributions to V2 in the t channel for two values of the scattering angle, 6 = 3° 
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with c = a'^/{TT'^s). The ^1,2,3 and Fa are de- 
fined as in [17 and [37j, and i?h(s) = 0{s — 

AMD i?had(s). ' 

Numerical results are given in Tables [3] and 
m where we include also the QED Born predic- 
tion, the Standard Model effective Born predic- 
tion for Mz = 91.188 GeV, Tz = 2.495 GeV, 



^^' = 0.23 and the contribution from the 
running of the fine-structure constant. The cut 
on the energy of the soft photons is set to a/s/2. 
We can see that, refering to the per mil accuracy: 
(i) electron vertices dominate over the rest of ver- 
tices (however, it is known that they largely can- 
cel with the contribution of the soft electron pair 
emission of [3], see also [31]); (h) contributions 
from infrared-safe boxes in Eq. [10] are substan- 
tial, mostly due to the factorizing diagrams; (iii) 
hadronic contributions play an important role. 

5. Summary 

Virtual NNLO QED corrections to massive 
Bhabha scattering have been completed in the 
small electron mass limit. Photonic, electron and 
heavy-fermion contributions have been checked 
by independent groups and different methods of 
calculations. Hadronic contributions have been 
calculated through the dispersion relation ap- 
proach; the kernels employed have been checked 
through a comparison with the heavy-fermion re- 
sult of HMS]. 
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Table 2 

Contributions to ReV2 in the s channel. See Ta- 
ble Q] for further details. 




The analysis of the virtual NNLO contri- 
butions shows that the results can influence 
Bhabha physics; therefore, further studies in- 
cluding Monte Carlo calculations with real 
bremsstrahlung are welcome [53]. 
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